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It is believed that quark matter can exist in neutron star interior if the baryon density 
is high enough. When there is a large isospin density, quark matter could be in a pion 
condensed phase. We compute neutrino emission from direct Urea processes in such 
a phase, particularly in the inhomogeneous Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) 
states. The neutrino emissivity and specific heat are obtained, from which the cooling 
rate is estimated. 



1. Introduction 

The phase structure of quantum chromodynamics (QCD) is one of the most chal- 
lenging problem in particle and nuclear physics. We schematically illustrate in Fig. 1 
the phase diagrams in T—fiB and T—p^i plots, where T, /ib and fj,j are the tempera- 
ture, baryon and isospin (or equivalently electron) chemical potentials respectively. 
In the T — hb diagram, the left panel of Fig. 1, the hadronic phase locates at low 
T and low ft b region and undergoes a phase transition or a crossover to the decon- 
fined quark phase at certain critical temperature T c or baryon chemical potential 
[iBc of the orders of T c ~ 200 MeV or fiBc ~ 1 GeV. At very high temperature, the 
quark-gluon-plasma (QGP), made of free quarks and gluons, forms. At asymptoti- 
cally high [is but low T, the ground state of QCD is the color- flavor-locked (CFL) 
superconductor 1 where the condensation of quark pairs spontaneously breaks color 
and chiral symmetries. At intermediate T and hb, although quarks and gluons are 
deconfined they are still strongly coupled. In this regime, many QCD phases are 
proposed in recent years, such as, at low T, two-flavor color superconductivity (2SC) 

2 , gapless 2SC (g2SC) 3 , gapless CFL (gCFL) 4 , spin-1 color superconductor 5 ' 6,7 , 
kaon condensation in the CFL phase 8 , et al.. For reviews of color superconductivity, 
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see, e.g. Ref. 9 . There may exist resonance states at intermediate T, such as strongly 
coupled QGP (sQGP) 10 at low \i B or the pseudo-gap phase at moderate \iB 11 • In 
the T—fii diagram, the right panel of Fig. 1, the hadron phase is in the region with 
low T and low fij while the QGP phase locates at very high T. At low T, when 
/zj grows above the value of the pion mass m v , the ground state turns out to be 
a Bosc-Einstein condensation (BEC) of pions, as fii increases further and is larger 
than about 230 MeV, the pion BEC crossover smoothly into the BCS superfluid of 
quark-anti-quark pairs with the condensate (uij^d) or (dij^u) 12 . At intermediate 
T, resonance states such as sQGP may occur at low fij while the states with strong 
fluctuations of thermally excited mesons or Cooper pairs are possible at moderate 

in. 




Fig. 1. The schematic phase diagrams of QCD on T — /j,b and T — /ij planes. 



In this paper we consider the quark matter cores in neutron stars. The neutrino 
emission from direct Urea processes d^>u + e~ + 9,u + e~ d + v is the most 
efficient way of cooling in quark matter. Schafer and Schwenzer summarized the 
neutrino emissivities and specific heats for a variety of color superconducting phases 
of quark matter 13 . Due to beta equilibrium the isospin chemical potential \ij is 
nonzero, quark matter could be a pion BEC (fii < 230 MeV) or a BCS superfluid 
(hi > 230 MeV) when /ij > 12 . On the other hand, the baryon chemical 
potential is still large (in the order of 1 GeV) and makes a big mismatch between 
the fermi surfaces of the pairing quarks u(u) and d(d), thus the BEC or BCS state 
is gapless. Such a gapless phase is stable in the BEC region but unstable in the BCS 
one with respect to the formation of nonzero LOFF momentum. In a previous work 
14 , we have studied the neutrino emissivity and cooling rate due to Urea processes 
for the gapless pion condensed quark matter in the BEC region. In this paper we 
study the neutrino emission in the LOFF phase. We work in two-flavor case in the 
moderate baryon density, where the role of strange quarks is not important. 

Our units are H = k B = c = 1 except particular specifications. As a convention, 
we denote a 4-momentum as K 11 = (fc ,k), and its 3-momentum magnitude as 
k = Ikl. 
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2. Quark Propagator 

Our starting point is the two flavor Nambu-Jona-Lasinio Lagrangian of QCD 

C = tpii^df, + ^70 - m )tp + g[{i>il>) 2 + (#f*7 5 V>) 2 ], (1) 

where -0 = (ii,<i) T is the quark fields, g is the coupling constant and f is the 
Pauli matrices. We have introduced the chemical potential matrix in flavor space, 
(i = diag(^ u , fj, d ) = (fj, + 5fj,,fj,- Sfi) = (/i S /3 + ^s/3 - m/2) with fj, B , m 

the baryon and isospin chemical potential respectively. We assume that the [3- 
cquilibrium is reached fj, d = Hu + He, which gives fij = —fj, e . In the chiral limit 
and without the chemical potentials, the Lagrangian (1) respects the symmetry 
Ub{1) <8> SUv(2) <8> SUa(2) corresponding to the baryon number, isospin vector 
and pseudovector conservation respectively. However, the presence of the chemical 
potentials explicitly break the isospin symmetry down to LV(1) with the conserved 
quantum number of r 3 , and chiral symmetry down to E/a(1) with the conserved 
quantum number of 175X3. By introducing the chiral condensate a = —2g(tpip) and 
pion condensates ir~ = Ae -2 = 4g(uij5d), ir + — Ae = Ag(dij5u), we arrive 
at the mean field Lagrangian 

£mf = ^ fi^d, + ^70 -m .M l5 \ _ <^ + A^ 
\ it «7 5 i^d^ + n d j a -mj Ag 

with the effective quark mass m = mo + a. There are mismatches in Fermi surfaces 
of anti-u and d quarks or anti-d and u quarks by the baryon chemical potential, 
hence pion condensates with nonzero total momentum or the LOFF states may 
be favored. The formation of the condensate a ~ (V>'0) breaks the Ua(1) chiral 
symmetry spontaneously with the Goldstone boson tto ~ ^7573^, and that of pion 
condensates A <~ (dij^u) ~ (uij 5 d) break the Uy(l) isospin symmetry sponta- 
neously. The translational and rotational symmetries are spontaneously broken by 
nonzero LOFF momentum 1. The partition function of the system can be written 
as a functional integral 



/ 




exp dr d 3 x£ mf . (3) 



Rewriting the quark fields via a gauge transformation which leaves the partition 
function unchanged, Xu{x) = u(x)e~ lhx , Xd(x) — d(x)e (we still call \u.d quark 
fields), the inverse quark propagator in flavor and momentum space reads 

V A *7s 7^ + 1-7 + ^70 -m) ' 

Note that K M represents the 4- momentum of the x fields instead of the ip fields. 
The propagator is written as 

( S UU (K) S ud (K)\ 
S(K)= {s du (K)S dd (K))- (5) 
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A straightforward calculation from Eq. (4) gives the four elements 



Suu(K) 
S dd (K) 
Sud(K) 

Sdu(K) 



{^K +ll + m)(K 2 _ - m 2 - A 2 ) + 2A 2 7 % 



(K 2 



m 2 + A 2 )(K 2 - m 2 + A 2 ) - A 2 [(K + + KJ) 2 
{^K^ + m){K 2 -m 2 - A 2 ) - 2A 2 7 % 



Am 2 }' 



(K 2 - m 2 + A 2 )(K 2 



A 2 ) - A 2 [(K+ + K_) 2 - Am 2 } ' 



(^K+^ + m)(K 2 _ -m 2 



A 2 



2A V? 



(K 2 



K 2 



A 2 ) (if 2 



A 2 ) - A 2 [(iT + + K-) 2 - Am 2 } 



{^K-^ + m){K 2 + - m 2 - A 2 ) - 2A 2 7 ^ 



(K 2 



K% 



A 2 ) (if 2 - m 2 + A 2 ) - A 2 [(if+ + K-) 2 - Am 2 } 

71 ■ <{ A 

-*7 A, 



(6) 

where if± = (fco+/i±<fyx, k±l) and l M = (Sfi, 1). The excitation spectra of the quasi- 
particles can be obtained by solving the equation det S~ 1 (ko, k) = or equivalently 
the roots of the denominator of the propagator for ko, 







'2 



[(**> + 
x[(fc 



m 



-m 2 + A 2 )(K 2 _- 
M + ^ ; ,) 2 -«, + ^) 2 -A 2 ] 
' 2 ^) 2 - A 2 



A 2 )-A 2 [(if+ + if_) 2 

2 A 21 



4m 2 



(7) 



where e k j = (i?k+i ± .Ek-i)/2 with i?k = Vk 2 + to 2 . To arrive at the last line, wc 
have taken the assumption that both A and 1 are small comparing to the quark 
Fermi momenta in the LOFF phase. We have four excitation branches, £> (k, 1) = 

— ry (e k { — a5/.i) 2 + A 2 — (p — as^j), with a, r = ±. Taking the same approximation 
to the numerators in the elements of the quark propagator, we neglect the terms 
proportional to A 2 1. Thus we rewrite the elements of the quark propagator as, 

B?(k,l)A£ +l7o 



S UU (K)~ ]T 



a.r— ± 



S M (K)~ J2 



a.r— ± 



k -E?(k,l) ' 

i?° r (k,l)A k ° l7o 
k -E?(k,l) ' 



where we have introduced the energy projectors, 



AO 

iv k 



1 + a 



7o ( 7 -k + m) 



and the Bogoliubov coefficients, 



S, a (k,l) = ^ 



£k 



k.l 



adfi 



(e+ - a^/i) 2 + A 2 



(8) 



(9) 



(10) 
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3. Neutrino Emissivity 



Since the characteristic energy scale of the Urea processes is much lower than the 
W-boson mass, we can use the Fermi current-current interaction to describe the 
Urea processes, 

Ant = ^=J"jt, (11) 

where the weak currents are defined by 

J*{x) = u-f(l - 75 )e + e-^Xul^l - 75)Xd, 

4( x ) = ^ - ^> + eM ' X X^(l - 75)X«, (12) 

where i and i for i = e,v denote spinor fields for electrons and neutrinos. Here 
G = GpCOsOc ~ 1.13488 x 10~ n MeV~ 2 is the four-fermion coupling constant. 

In the /3-equilibrated quark matter, the neutrino emissivity, defined as the total 
energy per unit time and per unit volume carried away by neutrinos and anti- 
neutrinos in escaping a neutron star can be written as 14 , 

xL X ^{P v ,P e )ImU^{E e -He-E v ,p e -p v ), (13) 

Here the on-shell 4-momenta for leptons is denoted by Pi — (E{,pi) (i — e,v), where 
the energies are Ei = E pi = y/pf + mf with m v — and m e w 0. The factor 2 in 
the front is due to the approximately identical cross sections for /3-decay and elec- 
tron capture processes at very low temperatures far from the condensation-normal 
transition temperature T c . But as the temperature approaches T c , one cannot take 
the contributions from the two processes as equal, as is demonstrated in Ref. 14 . We 
used Ub(x) = (e x l T — 1) _1 and Uf(x) = (e x l T + 1) _1 to denote the Bose-Einstein 
and Fermi-Dirac distribution functions, which maintain the detailed balance and 
Pauli enhancing/blocking effects for electrons in the electron capture//3-decay pro- 
cesses. Note that there are two identical terms in evaluating the emissivity (13), see 
Fig. 3 of Ref. 15 , which gives the correct overall factor and can easily be overlooked. 
We have dropped the chemical potential of neutrinos since there is no accumula- 
tions of neutrinos at typical temperatures inside an aging star. The leptonic tensor 
reads 

L x ° {P v , P e ) = Tr[ 7 A (l - 75 ) 7 • P e7 CT (l - 75 ) 7 • P v ] 

= 8[P e A P; + P A P; - P e • P v g x ° - ie Xa ^P ea P^}. (14) 
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The W-boson polarization tensor U Xcr (q ,q) can be written as 

d 3 p u n F [E?(p u , 1)] - n F [E- b (p d , 1)] 



a.b.r.s 



(2n)HE u E d qo - Es\p d , 1) + E?(p u , 1) 



xB r "(Pu, l)BZ a s (p d , 1)H%(P U , P d ), (15) 

where p d0 = p u o+qo, Pd =p„+q+21, P u = (E u ,ap u + al) and P d = (E d ,bp d -b\) 
with E u = E Pu+ \, E d = E Pd _i. Here p„, p d are the momenta of Xu, Xd- The bosonic 
and fermionic Matsubara frequencies are given by qo = i2mwT and p u o = i(2n + 
l)nT (m,n are integers). The quark tensor H x £ (P u , P d ) is defined by 

H%(P u ,P d ) = 4£„£ d Tr[7 A (l - 75 )A^ +1 7o7 CT (l - 7 5 )A^_ l7 o] 

= 8[P X P! + P X P^ - P u ■ P d g Xa - ie Xa ^P ua P d0 }. (16) 

By an analytic extension for the Matsubara frequency qo — i2mirT —* qo + i0 + , 
where the second qo is real, the imaginary part of the retarded polarization tensor 
of W-bosons can be read out directly, 

d 3 p u 



Imll^o, q) = ^ E / (2n)HE u E d S[q ° ~ 1} + KiPu ' Pd) 

B a r {p u ,\)BZ b s {p d ,\). (17) 



a.b,r,s * 

n F [E?(p u ,l))n F [-E- b (p d ,l)] oa ,_ 



n B (-qo) 

Substituting Eq.(17) into Eq.(13) we arrive at 



4iV n ^ 



a.b.r.s 



d 3 p e d 3 p u d 3 p u d 3 p d 



(2w) 3 2E e (2w) 3 2E„ {2ir) 3 2E u {2nf2E d 

xE u {2^f5[E e -ii e -E v + E a r {p u , 1) - E~ b (p d , 1)] 

x<5 3 (p e - p v + p u - p d + 2\)B a r {p u , l)BI b s (p d , 1) 

xn F (E e - Ve)n F [E?(p u ,l)]n F [-E7 b (p d ,l)]\M ab \ 2 . (18) 

We have introduced the shorthand notation |M a fc| 2 , the spin-averaged scattering 
matrix element of /3-decay or electron capture 19 , 

\M ab \ 2 = ^-L Xa (P e , P u )H x Z(P u ,P d ) = 64G 2 (P e • P U ){P V ■ P d ), (19) 

where the energy projection indices a, b are hidden in the quark momenta. 

In the LOFF phase, chiral symmetry is almost restored, so we can safely set 
quark masses zero, m Uyd — 0. We consider the low temperature and high chemical 
potentials for quarks, the dominant contribution of the phase space integral in quark 
momenta comes from the gapless modes of the excitations with positive energies, 
i.e. terms with r = s = — . Also, since the gap parameter and the LOFF momentum 
are small comparing to the quark chemical potentials, the Bogoliubov coefficients 
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with a = — and b = — are strongly suppressed, so we only keep the term with 
a = + and b = +. For electrons, the relevant momenta are near the chemical 
potential. Recalling that E+(p u ,l) and £JZ(pd,l) are the dispersion relations for 
quasi- u-quarks and quasi-d-quarks, one can obtain the gapless momenta as 

P°u/ d »(1-k)[^t1 cos6 u/d ) 2 - A 2 ± Sfi] (20) 

for u/d quarks, where k is due to the Fermi-liquid correction. In perturbative QCD 
we have k = 2a s / (3n) with a s being the strong coupling constant 17 . In NJL model 
we have 18 n = 4g \j? B / (3tt 2 ) . We have set the z-direction along the LOFF momentum 
1 and denote as 0i and the polar and azimuthal angles of pj, i = e,v,u,d. Near the 
gapless momenta, we can make expansion E±(p u , 1) w v u \p u — p^\ and EZ(pd, 1) ~ 
Vd\Pd — P d \ with velocities v u / d = ^/l — A 2 / (Ji =p Z cos 9 u / d ) 2 . The matrix element is 
now evaluated as 

|M ++ | 2 w MG 2 p ePv E u E d [l - (1 - K )cos0„ e ][l - (1 - k)cos^], (21) 

where fl^f and ue are angles between p d — 1 and p„ and between p u + 1 and p e . 
Since both A, 1 are small and p v <~ T -C n, 5/j, is negligible, the delta function for 
the energies in Eq. (18) can be rewritten in the form 6[E e — fi e — E v + EZ(p u , 1) — 
El(p d ,l)} w p e /(p°p a d )5(cos9 ud - cos#° d ) with 6 ud the angle between p d - 1 and 
p u + 1 and cos0O d w [(p°) 2 + {p°) 2 - 2/j, e l(cos6 u + cos6 d )}/(2p u p° d ). The integration 
over electron three-momentum can be carried out by the delta function about the 
momentum conservation. Now we get, 

5" NcG MeT 7 (2tt)3 (2tt)3 (2 W )3^ B - (P « )B + (P " ) 
x<5(cos#„ d - cos6»° d )[l - (1 - k) cos#„ e ][l - (1 - k) cosd dv \ 

= ^N C G 2 ^ f ^^MsttfjB^ 
315 J (2ir) 6 (2it) 6 v u v d 

x <5(cos 9 ud - cos 6° ud ) [1 - (1 - k) cos 9 ue ] , (22) 

where we have used 16 f °° dp v dp u dp d p%n F (p v + \Pd - P° d \ - \Pu - p"|)«f(|p« - 
P°u\) n F(-\Pd ~P d \) « 457tt 6 T 6 /5040. By setting A = I = in Eq. (22), one can 
recover the well-known result for the neutrino emissivity in normal quark matter 

19 

e " M 7^* N cG 2 ^u»d ( 1 + — ) ^T 6 . (23) 
2520 \ /U„/ 

We make numerical evaluation of the neutrino emissivity. We choose k = 2/{3n), 
li = 300 McV, —<5/i = 125 MeV. These are typical values to support the possible 
LOFF pion superfluid. In Fig. 2, the neutrino emissivity in unit of eo is presented as 
function of the pion condensate A and the LOFF momentum I. Due to the gapless 
nature, along the A axis, these is no exponential but an approximate quadratic 
suppression. The appearance of the LOFF momentum also slightly lowers the neu- 
trino emissivity. Note that the approximation we made here is somewhat different 



914^ 
315'' 



22, 2009 10:10 WSPC/INSTRUCTION FILE LOFF'UrcaTion 



8 Huang, et.al. 




Fig. 2. (Color online) The neutrino emissivity in a LOFF pion superfluid as functions of the pion 
condensate A and the LOFF momentum I. The parameters are set to k = 2/(3n), fi = 300 MeV, 
-Sfj, = 125 MeV. 



from that we did in the previous work 1 . In this paper we evaluate the Bogoliubov 
coefficients at gapless momenta and result in a suppression in the emissivity. In the 
previous work the Bogoliubov coefficients are taken to be 1 approximately which 
leads in turn to an emissivity enhancement. The latter approximation is valid even 
at high temperatures, but the former one should be adopted only at low tempera- 
tures. 



4. Cooling Rates 

To get the cooling behavior, one must know the specific heat of the pion superfluid 
which can be derived from its definition cy(T) = T(dS / dT)v with S being the 
entropy density, 

^{n F [E?(k,l)]lnn F [E?(k,l)] 

r,a 

+n F [-E a r (k, 1)] In n F [-E* (k, 1)] } , (24) 

where the pre-factor 2N C comes from the degeneracies of spins and colors. Since 
we work at low temperatures and high quark chemical potentials, the dominant 
contribution to the specific heat are from the gapless momenta of positive energy 
excitations, 



cy 



(T)«2iV c ]T J ^nr[ElQL,l)]n F [-E°QL,l)] 



(2tt) 3 ri - y ' n r 1 — ^'-' J T2 
2k 2 N c ^ f cM (rf) 2 
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where we have used the fact that at low temperatures A, 1 are almost constants. 
When A and / vanish, we reproduce the specific heat of normal two-flavor quark 
matter c vo = 7T with 7 = 2N c (pl + fJ? d )T/6. 

The time evolution of the temperature can be obtained by solving the following 
equation t - t = - dT' c v (T') / 'e(T') , where To is the temperature at an initial 
time to- Substituting Eq. (22) and Eq. (25) into this equation, we arrive at 

rT dT'G(AJ) 



t-t 



210 



4577r 3 G 2 /z, 
where we have introduced the notations 



/ 



T' 5 F(A,0' 



(26) 



dfl u dfl d p° uPd 



F(A,l) = 

•—u.d 



(2tt) 3 (2tt)3 v u v d 
dfl (p°) 2 



B±( P o u )B-( P ° d )6(cos0 ud - cosOtl - (1 - «)cos0 Me ] 

(27) 




10 20 30 40 50 

t(yr) 

Fig. 3. The cooling curves of the pion supcrfluid in the LOFF phase and of normal quark matter 
with To = 0.1 MeV at to = 1 yr. The parameters are the same as in Fig. 2. 



Fig. 3 shows the time evolution of the pion supcrfluid in the LOFF phase with 
the initial temperature T = 0.1 MeV at to = 1 yr. For comparison, we also show the 
curve of normal quark matter. We choose A = I — 50 MeV. The other parameters 
are chosen as the same as in Fig. 2. We see that the cooling of the LOFF pion 
supcrfluid is slightly slower than that for normal quark matter. This situation is 
very similar to the case of inhomogeneous CFL quark matter, see 20 . 

5. Summary 

If the isospin (or equivalently electron) chemical potential in a neutron star is large 
enough, the LOFF pion superfiuid can be a possible ground state in the interior 
matter of the neutron star. We calculated the neutrino emissivity, specific heat and 
then the cooling rate through direct Urea processes in this phase. The main results 
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are shown in Eq. (22) for the neutrino emissivity and in Eq. (26) for the cooling 
rate. At low temperatures, the neutrino emissivity is slightly smaller in the LOFF 
pion superfluid than that in normal quark matter, as indicated in Fig. 2 which leads 
to a slightly slower cooling rate for the LOFF phase as shown in Fig. 3. 
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It is believed that quark matter can exist in neutron star interior if the baryon number 
density is high enough. When there is a large isospin density, quark matter could be 
in a pion condensed phase. We compute neutrino emission from direct Urea processes 
in such a phase, particularly in the inhomogeneous Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) states. The neutrino emissivity and specific heat are obtained, from which the 
cooling rate is estimated. 

1. Introduction 

The phase structure of quantum chromodynamics (QCD) is one of the most chal- 
lenging problem in particle and nuclear physics. We schematically illustrate in Fig. 
1 the phase diagrams in T — fis and T — fij plots, where T, ^b and are the 
temperature, baryon and isospin (or cquivalently electron) chemical potentials re- 
spectively. In the T—fis diagram, the left panel of Fig. 1, the hadronic phase locates 
at low T and low [Ib region and undergoes a phase transition or a crossover to the 
deconfined quark phase at certain critical temperature T c or baryon chemical poten- 
tial hbc of orders of T c ~ 200 MeV or hb c ~ 1 GeV. At very high temperature, the 
quark-gluon-plasma (QGP), made of free quarks and gluons, forms. At asymptoti- 
cally high [is but low T, the ground state of QCD is the color- flavor-locked (CFL) 
superconductor 1 where the condensation of quark pairs spontaneously breaks color 
and chiral symmetries. At intermediate T and hb, although quarks and gluons are 
deconfined they are still strongly coupled. In this regime, many QCD phases are 
proposed in recent years: two- flavor color superconductivity (2SC) 2 , gapless 2SC 
(g2SC) 3 , gapless CFL (gCFL) 4 , spin-1 color superconductor 5 ' 6 ' 7 , kaon conden- 
sation in the CFL phase 8 , et al.. For reviews of color superconductivity, see, e.g. 

* e-mail: huangxg03@mails.tsinghua.cdu.cn 

t e-mail: qunwang@ustc.edu.cn 

t e-mail: zhuangpf@mail.tsinghua.edu.cn 
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Ref. 9 . There may exist resonance states at intermediate T, such as strongly cou- 
pled QGP (sQGP) 10 at low /is or the pseudo-gap phase at moderate hb 11 • In 
the T — \ii diagram, the right panel of Fig. 1, the hadron phase is in the region 
with low T and low fii while the QGP phase locates at very high T. When fii 
grows above the value of the pion mass m ff , the ground state turns out to be a 
Bose-Einstein condensation (BEC) of pions, as /xj increases further and is larger 
than about 230 MeV, the pion BEC crossover smoothly into the BCS superfluid of 
quark-anti-quark pairs with the condensate (ui^d) or (di^u) 12 . At intermediate 
T, resonance states such as sQGP may occur at low /ij while the states with strong 
fluctuations of thermally excited mesons or Cooper pairs are possible at moderate 
Hi. 




Fig. 1. The schematic phase diagrams of QCD on T — /j,b and T — pi planes. 



In this paper we consider the quark matter cores in neutron stars. The neutrino 
emission from direct Urea processes d — > u + e~ + P,u + e~ — > d + is the most 
efficient way of cooling in quark matter. Schaefer and Schwenzer summarized the 
neutrino emissivities and specific heats for a variety of color superconducting phases 
of quark matter 13 . In two- flavor quark matter, due to beta equilibrium the isospin 
chemical potential fii is nonzero, quark matter could be a pion BEC (/ij < 230 
MeV) or a BCS superfluid (pi > 230 MeV) when fii > m„ 12 . On the other hand, 
the baryon chemical potential is still large (in the order of 1 GeV) and makes a big 
mismatch between the fermi surfaces of the pairing quarks u(u) and d(d), thus the 
BEC or BCS state is gapless. Such a gapless phase is stable in the BEC region but 
unstable in the BCS one with respect to the formation of nonzero LOFF momentum. 
In a previous work 14 , we have studied the neutrino emissivity and cooling rate due 
to Urea processes for the gapless pion condensed quark matter. In this paper we 
study the neutrino emission in the LOFF phase. We work in two-flavor case in the 
moderate baryon density, where the role of strange quarks is not important. 

Our units are h = ks = c = 1 except particular specifications. As a convention, 
we denote a 4-momentum as = (fco,k), and its 3-momentum magnitude as 
k = |k|. 
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2. Quark Propagator 

Our starting point is the two flavor Nambu-Jona-Lasinio Lagrangian of QCD 

C = tp(i^dft + ^70 - m )V> + g[(^) 2 + (#ry 5 V>) 2 ], (1) 

where = {u,d) T is the quark fields, g is the coupling constant and f is the 
Pauli matrices. We have introduced the chemical potential matrix in flavor space, 
/x = diag(^ u , fj, d ) = (fj, + 6fi,fi- 6fi) = (/xb/3 + /J-i/2, ^ B /3 - M//2) with /zj 
the baryon and isospin chemical potential respectively. We assume that the /?- 
cquilibrium is reached fid — fi u + /j, e , which gives fii = —fi e . In the chiral limit 
and without the chemical potentials, the Lagrangian (1) respects the symmetry 
C^s(l) <8> S£V(2) <g> SUa(2) corresponding to the baryon number, isospin vector 
and pseudovector conservation respectively. However, the presence of the chemical 
potentials explicitly break the isospin symmetry down to fv(l) with the conserved 
quantum number of T3, and chiral symmetry down to Ua{{) with the conserved 
quantum number of 175X3. By introducing the chiral condensate a — —2g(tptp) and 
pion condensates ir~ — Ae _2llx = Ag{ui^^d), ir + = Ae 2 * lx = 4g(di^u), we arrive 
at the mean field Lagrangian 

r 7 ( il^d^ + fiujo ~ m TT+i"/ 5 \ a 2 + A 2 



7r «7 5 «7^ + Md7o - m J Ag 

with the effective quark mass m = m + a. There are mismatches in Fermi surfaces 
of anti-u and d quarks or anti-d and u quarks by the baryon chemical potential, 
hence pion condensates with nonzero total momentum or the LOFF states may 
be favored. The formation of the condensate a <~ (V>'0) breaks the J7a(1) chiral 
symmetry spontaneously with the Goldstone boson ttq <~ tpi^Tstp , and that of 
pion condensates A ~ (di^u) ~ (ui^d) break the f/y(l) isospin symmetry spon- 
taneously with the Goldstone boson ao ~ i>T3i[). The translational and rotational 
symmetries are spontaneously broken by nonzero LOFF momentum 1. The partition 
function of the system can be written as a functional integral, 



J [#][#]< 




cxp dr d 3 x£ mf , (3) 



Rewriting the quark fields via a gauge transformation which leaves the partition 
function unchanged, Xu( x ) = u(x)e~ ti ' x , Xd{x) — d(x)e d x (we still call \u,d quark 
fields), the inverse quark propagator in flavor and momentum space reads 

s -i (K \ = / 7^^ - 1 • 7 + M«7o - m Ai 75 \ 

1 ' V A *7s YK„ + 1 • 7 + (idlo -mj' [ ' 

Note that represents the 4- momentum of the \ fields instead of the ip fields. 
The propagator is written as 
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A straightforward calculation from Eq. (4) gives the four elements 

{^K +ll + m)(K 2 _ - m 2 - A 2 ) + 2A 2 7 % 



Suu(K) 
S dd {K) 
Sud(K) 

S du (K) 



{K% -m 2 + A 2 ){K 2 - m 2 + A 2 ) - A 2 [(K+ + KJ) 2 - 4m 2 ] ' 

{^K^ + m){K\ - m 2 - A 2 ) - 2A 2 7 % 
(K% -m 2 + A 2 )(K 2 -m 2 + A 2 ) - A 2 [(K + + KJ) 2 - Am 2 ] ' 

{^K +ll + m)(K 2 - m 2 - A 2 ) + 2A 2 7 ^ 

(K\ -m 2 + A 2 )(K 2 -m 2 + A 2 ) - A 2 [(K+ + K_) 2 - Am 2 ] 

Kt - m 2 

{^K-^ + m){K\ - m 2 - A 2 ) - 2A 2 7 % 
(Kl -m 2 + A 2 )(K 2 _ -m 2 + A 2 ) - A 2 [(K+ + KJ) 2 - Am 2 ] 

x 7-^2 9 n 5 A, 



(6) 

where K± = (fco + /i±<5/i, k±l) and — (6/j,, 1). The excitation spectra of the quasi- 
particles can be obtained by solving the equation det S~ 1 (ko, k) = or equivalently 
the roots of the denominator of the propagator for ho, 



Q=(Kl-m 2 + A 2 )(K 2 _ - m 2 + A 2 ) - A 2 [(K+ + KJ) 2 - Am 2 ] 
« [(fco+^ + £ k ,,) 2 -(<, + 5/i) 2 -A 2 ] 

x[(fc +M-£k,i) 2 -«i-^) 2 -A 2 ], (7) 



where = (-Ek+i ± £ ; k-i)/2 with = y/k 2 + m 2 . To arrive at the last line, we 
have taken the assumption that both A and 1 are small comparing to the quark 
Fermi momenta in the LOFF phase. We have four excitation branches, 



£ r a (k, 1) = -ry - aS^) 2 + A 2 - (/x - ae k l ), (8) 

with a,r = ±. Taking the same approximation to the numerators in the elements 
of the quark propagator, we neglect the terms proportional to A 2 1. Thus we rewrite 
the elements of the quark propagator as, 

M ] ~ a ^ ± fco-^(k,l) ' 

5° r (k,l)A k ° l7o 
^ fc -^(k,l) ' (9) 

a,r=± ' v ' 

where we have introduced the energy projectors, 

7o ( 7 -k + m) 



A? = I 



1 + a- 



(10) 
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and the Bogoliubov coefficients, 



5 r a (k,l) = - 



1 — ar 



-k,i 



aSfj, 



^(e+ 1 -a^)2 + A2_ 



(11) 



Note that Sfj, = (fi u — /Ud)/2 here, which differs by a sign from that in Ref. 14 , so 
the Bogliubov coefficients in Eq. (11) have the same convention in signs of a, r as in 
Eq. (12) of Ref. 14 . We see that the propagator elements in Eq. (9) are in a similar 
form to Eq. (8,9) of Ref. 14 . There are some differences: (i) The energy excitation 
spectra in the current case are slightly different from that in Ref. 14 by a shift ae_ 
and by that the energy is e£ Y instead of E^; (ii) The Bogliubov coefficients depend 
on e£ Y not on E^; (iii) The energy projectors for u and d quarks depend on k ± 1 
instead of both on k. 



3. Neutrino Emissivity 

Since the characteristic energy scale of the Urea processes is much lower than the 
W-boson mass, we can use the Fermi current-current interaction to describe the 
Urea processes, 

Ant = j=J"Jt, (12) 

where the weak currents are defined by 

J*(x) - - 75 )e + e-^Xul^l - 75)Xd, 

= g ^(l ~ 7*)" + e M ' x x^(l - 75)X«, (13) 

where i and i for i = e, v denote spinor fields for electrons and neutrinos. Here 
G = GfCOsOc ~ 1.13488 x 10 _11 MeV~ 2 is the four-fermion coupling constant. 

In the /3-equilibrated quark matter, the neutrino emissivity, defined as the total 
energy per unit time and per unit volume carried away by neutrinos and anti- 
neutrinos in escaping a compact star. Using the Lagrangian (12) in the Kadanoff- 
Byam equation, the neutrino emissivity can be derived as 14 , 

xL Xa (P u ,P e )ImUt(E e -n e -E v ,p e -p v ), (14) 

Here the on-shell 4-momenta for leptons is denoted by Pi = (Ei,pi) (i — e,v), where 
the energies are Ei = E Pi = y/pj + rrif with ra v = and m e w 0. The factor 2 in 
the front is due to the approximately identical cross sections for /3-decay and elec- 
tron capture processes at very low temperatures far from the condensation-normal 
transition temperature T c . But as the temperature approaches T c , one cannot take 
the contributions from the two processes as equal, as is demonstrated in Ref. 14 . We 
used tib(x) — (e x ' T — 1) _1 and Uf{x) — (e x / T + 1) _1 to denote the Bose-Einstein 
and Fermi-Dirac distribution functions, which maintain the detailed balance and 
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Pauli enhancing/blocking effects for electrons in the electron capture//3-decay pro- 
cesses. Note that there are two identical terms in evaluating the emissivity (14), see 
Fig. 3 of Ref. 15 , which gives the correct overall factor and can easily be overlooked. 
We have dropped the chemical potential of neutrinos since there is no accumula- 
tions of neutrinos at typical temperatures inside an aging star. The leptonic tensor 
L Xa (P e ,P v ) reads 

L x ° {P v , P e ) = Tr[ 7 A (l - 75 ) 7 • Pel a {l ~ 75)7 ' Pu] 

= 8[P A P; + P x P a e - P e • P v g x ° - ie Xa °Pp ea P v0 \. (15) 

The W-boson polarization tensor H Xa {q ,c\) can be written as 

n x ^ qo , q )^N c Tj2 1 ^|Tr[ 7 A (l- 75 )^ u (P ll ) 7 -(l- 75 )5 dd (P d )] 

^ivrVV f d3pu B ^ Pu > l) B -sM 

C naArj i^f PuO ~ E? ( Pu ,l) p M - E^ ( Pd ,l) 

xTr[ 7 A (l - 7 5)A^ +l7o7 CT (l - 75)A^_i7o] 

d 3 p u n F [P r a (p„, 1)] - n F [E- b (p d , 1)] 



= Nc j: 



,b,r,; 



(2n)HE u E d qo - p- b (p d , 1) + E?(p u , 1) 



xP r a (p M , \)BZ b s (p d , \)H%{P U , P d ), (16) 

where P u = (p„ ,P«), P d = (Pdo,Pd) and Q = (q ,q) which satisfy P d = P u + 
Q + (0,21). Here p«,p,i are the momenta of Xu,Xd- The bosonic and fermionic 
Matsubara frequencies are given by q n = i2rmtT and p u0 = i(2n + 1)ttT (m, n are 
integers). The retarded polarization tensor II ^ is obtained from the above after 
taking the analytic extension for the Matsubara frequency qo = i2mirT — > qo +i0 + , 
where the second qo is real. The quark tensor H X £(P U , P d ) is defined by 

H x £(P u ,P d ) = 4P u P d Tr[ 7 A (l - 75)A^ +1 7o7 ff (l - 7 5 )A^_ l7 o] 

= 8[P A P d CT + P X P: - P u ■ P d g x ° - ie Xa ^P ua P d0 ] , (17) 

where P u = (E u , ap u + al), P d = (E d , bp d -b\) with E u = E Pu+h E d = P Pd -i. Using 
ImII A i cr (<7o, q) = ImiT A<T (<7o +«0 + , q), the imaginary part of the retarded polarization 
tensor of W-bosons can be read out directly, 

ImII*f fa,, q) =nN c £ J ^^L^- 5 [q E;»(p d , 1) + E?(p u , l)]H% (P u , P d ) 

a.b.r,s 

n F [E?(p u , l)]n F [-E- b (p d , 1)] Uft -b (n n 

x ( t B r (p u ,l)B_ s (p d ,l). (18) 
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Substituting Eq.(18) into Eq.(14) we arrive at 

d 3 p e d 3 p v d 3 p u d 3 p d 



4iV c £ / 

t, „ „ J 



{2ir) 3 2E e {2tt) 3 2E v {2ir) 3 2E u {2i:) 3 2E d 



a,b,r.s 1 

xE v (2n)*5[E e -^ e -E v + E a r {p u , 1) - E7 b (p d , 1)] 
x<5 3 (p e - p v + p u - p d + 21)B r a (p„, l)BZ b s (p d , 1) 
n F {E e -n e )n F [E^{p u ,\)]n F [-E-\p d ,\)]\M ab \ 2 . (19) 

We have introduced the shorthand notation |M a fc| 2 , the spin-averaged scattering 
matrix element of /3-decay or electron capture 19 , 

\M ab \ 2 = ^L Aa (P e , P v )H^{P u ,P d ) = 64G 2 (P e • P U ){P V ■ P d ), (20) 

where the energy projection indices a, b are hidden in the quark momenta. 

In the LOFF phase, chiral symmetry is almost restored, so we can safely set 
quark masses zero, m Utd — 0. We consider the low temperature and high chemical 
potentials for quarks, the dominant contribution of the phase space integral in quark 
momenta comes from the gapless modes of the excitations with positive energies, 
i.e. terms with r = s = — . Also, since the gap parameter and the LOFF momentum 
are small comparing to the quark chemical potentials, the Bogoliubov coefficients 
with a = — and b = — are strongly suppressed, so we only keep the term with 
a = + and b = +. For electrons, the relevant momenta are near the chemical 
potential. Recalling that i£+(p u ,l) and EZ(p d ,l) are the dispersion relations for 
quasi-u-quarks and quasi-d-quarks, one can obtain the gapless momenta as 

P°u/ d « (1 - «)[^ T I cos 6 u/d y - A 2 ± SfM] (21) 

for u/d quarks, where k is due to the Fermi-liquid correction. In perturbative QCD 
we have k = 2a s /(37r) with a s being the strong coupling constant 17 . In NJL model 
we have 18 k — 4g fi 2 B / \3tt 2 ) ■ We have set the z-direction along the LOFF momentum 
1 and denote as 9i and <f>i the polar and azimuthal angles of pj, i = e,v,u,d. Near the 
gapless momenta, we can make expansion E^(p u ,l) ss v u \p u — p°| and EZ{p d , 1) ~ 
Vd\Pd~ P d \ with velocities v u / d = yjl — A 2 /(^ =p lcos9 u / d ) 2 . The matrix element is 
now evaluated as 

|M ++ | 2 « 64G 2 p e p„E u E d (l P°+ lcose " cos6ue){1 Pd' loo S 9 d 

&u &d 



&u &d 

' QAG 2 p e p v E u E d [l - (1 -k) cos<9 ue ][l -(!-«) cos9 dv ], 



(22) 



where 9 dv and 9 ue are angles between p d — 1 and p v and between p u + 1 and p e , or 
cos g = ( p " + l l ' Pe = ( pd ~ 1 cos e J?- ~(Pu + l cos 6> ») 2 ~ Pi (23) 

|p u + l|p e 2(p u + I COS 9 u )p e 
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Since both A, 1 are small and p v ~ T -C p,5p, is negligible, the delta function for 
the energies in Eq. (19) can be rewritten in the form 

5[E e - Me - + E±{p u , 1) - EZ ( Pd , 1)] 
w S[E e - n e - E v + v u \p u -p°\ - v d \p d -p d \] 
-» S[E e -fi e -E v + ( Pu - p° u ) - ( Pd - p° d )} 
= Pe/(p u P d )5(cos6 ud -cos6 ud ) (24) 

The second line is due to the fact that the integral is an even function of (p u /d—p u / d ) 
so that we can restrict the integration range to p u /d~P u / d > an d double the result. 
We have used v u , v d w 1 and 

cog q d _ (Pu + I cos e u ) 2 + {p d -l cos 9 d ) 2 - pi 

2(p u + 1 cos e u ) (p d - 1 cos e d ) 

c °S#°d = ( cos ^«d)|p u =pO, Pd =pO )J , e=Me . (25) 

The integration over electron three-momentum can be carried out by the delta 
function about the momentum conservation. Now we get, 

£ ~ 315" ^ MeT 7 (27r)3 (2tt)3 (^)3PuP^_(p u )B + (p d ) 
X(5(cos6» Md - cos6>° d )[l - (1 - k) cos6» ue ][l - (1 - k) cos0 di ,] 

xd(cose ud - cos6»° d )[l - (1 - K)cos6i„e], (26) 
where we have used 16 , 

/ dp v dp u dp d pln F {p v + \Pd-p° d \- \Pu -P u \) 
Jo 

xn F (\ Pu - p° u \)n F (-\p d - p° d \) « 457^ 6 T 6 /5040. (27) 

By setting A = I = in Eq. (26), one can recover the well-known result for the 
neutrino emissivity in normal quark matter 19 , 

£ ° ~ ^Jj 7rN - G2 ^^d (l + j^j nT 6 . (28) 

We make numerical evaluation of the neutrino emissivity. We choose k = 2/{3tt), 
p = 300 McV, —<5/i = 125 MeV. These are typical values to support the possible 
LOFF pion superfluid. In Fig. 2, the neutrino emissivity in unit of eo is presented as 
function of the pion condensate A and the LOFF momentum I. Due to the gapless 
nature, along the A axis, these is no exponential but an approximate quadratic 
suppression. The appearance of the LOFF momentum also slightly lowers the neu- 
trino emissivity. Note that the approximation we made here is somewhat different 
from that we did in the previous work 14 . In this paper we evaluate the Bogoliubov 
coefficients at gapless momenta and result in a suppression in the emissivity. In the 
previous work the Bogoliubov coefficients are taken to be 1 approximately which 
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Fig. 2. (Color online) The neutrino emissivity in a LOFF pion superfluid as functions of the pion 
condensate A and the LOFF momentum I. The parameters are set to k = 2/(3tt), fi = 300 MeV, 
-Sfj, = 125 MeV. 



leads in turn to an emissivity enhancement. The latter approximation is valid even 
at high temperatures, but the former one should be adopted only at low tempera- 
tures. 



4. Cooling Rates 

To get the cooling behavior, one must know the specific heat of the pion superfluid 
which can be derived from its definition cy(T) = T{dS / dT)y with S being the 
entropy density, 

S = - 2N cJ2 J p{'>F[£;(k,l)]ln nf K(k,l)] 

+n F [-E a r (k, 1)] In n F (k, 1)] } , (29) 

where the pre-factor 2N C comes from the degeneracies of spins and colors. Since 
we work at low temperatures and high quark chemical potentials, the dominant 
contribution to the specific heat are from the gapless momenta of positive energy 
excitations, 

Cy (T)«2A^ J J0_ np [E a (k; 1)]njr { _ E a (k; 1}] I)] 2 

~ 2N < E / 7^3 (Pi) 2 dkn F ( Vi \k - p°\)n F (- Vi \k - p^ f^'/^ 
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where we have used the fact that at low temperatures A, 1 are almost constants. 
When A and / vanish, we reproduce the specific heat of normal two-flavor quark 
matter c vo = 7T with 7 = 2N c {pl + ^ 2 d )T/6. 

The time evolution of the temperature can be obtained by solving the following 
equation, 

'-"-ly-^r- 

where To is the temperature at an initial time to- Substituting Eq. (26) and Eq. 
(30) into the above equation, we arrive at 

,_, 210 f T dV G(A,Q 

*°- 457^GV7to T*F(A,iy [62) 
where we have introduced the notations 



F(A,l) = J (2^3 ^^ Bl(pl)B+(p d )S(co S e ud - co S e° ud )[l - (1 - k) cosfl, 
(2tt)3 Vi 



i—u,d 




Fig. 3 shows the time evolution of the pion superfiuid in the LOFF phase with 
the initial temperature To = 0.1 MeV at to — 1 vr - F° r comparison, we also show the 
curve of normal quark matter. We choose A = I = 50 MeV. The other parameters 
are chosen as the same as in Fig. 2. We see that the cooling of the LOFF pion 
superfiuid is slightly slower than that for normal quark matter. 

5. Summary 

If the electron chemical potential in a neutron star is large enough (> 230 MeV), the 
LOFF pion superfiuid can be a possible ground state in the interior matter of the 



18, 2007 10:59 WSPC/INSTRUCTION FILE LOFF - UrcaTion03 



Neutrino Emission From Pion Condensed Quark Matter 11 

star. We calculated the neutrino emissivity, specific heat and then the cooling rate 
through direct Urea processes in this phase. The main results are shown in Eq. (26) 
for the neutrino emissivity and in Eq. (32) for the cooling rate. At low temperatures, 
the neutrino emissivity is slightly smaller in the LOFF pion superfluid than that in 
normal quark matter, as indicated in Fig. 2 which leads to a slightly slower cooling 
rate for the LOFF phase as shown in Fig. 3. 
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